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ABSTRACT 

by 

Mario A. Baltazar 



Rational design methods based on a statistical 
description of both the ship load and strength have been 
proposed by several authors in lieu of the traditional 
bending moment calculations. In general, the effect of 
horizontal moments on the primary strength has been 
neglected. It is shown in this paper, as several others 
have found, that the horizontal moments can be very 
significant. Two methods are described to obtain the 
magnitude of the equivalent combined moment. 

An example application is shown for a large tanker 
by use of the MIT seakeeping program. Wave moment trends 
with speed, sea state and wave heading are noted. Measured 
and predicted values of the vertical bending stress are 
compared. Finally, the parameters of the weibull distrib- 
ution to describe long-term wave moment amplitudes are 
derived. 
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= ship's length 
= wave height 
= significant wave height 
= vertical bending moment 
= horizontal bending moment 
= equivalent combined bending moment 
= vertical section modulus 
= horizontal section modulus 
= section modulus ratio 

= mean square value of the amplitude of the response 
- parameters of the Weibull distribution 
= mean value of the response amplitude 
= amplitude of the vertical moment component 
= amplitude of the horizontal moment component 
= ship's speed 
= effective wave length 
= response amplitude operator 
= root mean square 
= circular frequency 

= phase angle of the horizontal moment with respect to a 
wave with its crest amidship 

= phase angle of the vertical moment with respect to 
a wave with its crest amidship 

= phase angle of the equivalent combined moment 

= correlation coefficient of the vertical and horizonral 
bending moment 

vi 



= total or combined stress 

= vertical bending stress 

a„ = horizontal bending stress 
H 

X = wave length 

S„(w) = spectral density of the response 

I\ 

H(to) = RAO of the ship 
S(w) = sea spectral density 
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Chapter I 



INTRODUCTION 

The traditional method for determining the load for 
primary strength calculations consists of poising the ship 
on a trochoidal wave of length equal to that of the ship and 
wave height equal to one- twentieth (1/20) of the length. 

The static bending moment is then determined by the load- 
shear-moment relations from beam theory. There have been 
various proposals to modify the wave height (e.g., H = 0.6L^*^ 
or H = l.lv^) partly to correct for dynamic effects and partly 
because it is thought than an L/20 wave is too severe for the 
very large ships. Notwithstanding, this empirically derived 
procedure was considered entirely satisfactory until relatively 
recently and is still in prevalent use inasmuch as it gives 
conservative estimates. 

The appearance of large ocean-going tankers and bulk 
carriers and novel ship types gave impetus towards the search 
for a "rational" procedure for ship structural design. It 
was realized that the application of the standard bending 
moment calculations for these ships lead to unrealistic re- 
sults. At the same time new knowledge relating to ship 
structural design had accumulated and experimental and comp- 
utational facilities to advance this knowledge have improved 
tremendously . 
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Many marine-oriented societies and organizations 
undertook basic research to improve the accuracy of pre- 
dictions of sea loads and motions, and the development of a 
rational design method applicable for all ship types. Full 
scale and model tests data were collected. The random 
nature of both the loads on a ship and its capacity to re- 
sist these loads was recognized, so that probability theory 
was utilized to develop statistical models for the design of 
ships [1-14]*. The advantages of a probabilistic approach 
over the traditional deterministic calculations are discussed 
in reference C7]. 

In the conventional method, the primary strength of 
a ship is based solely on the longitudinal stresses caused 
by vertical bending, or bending about the horizontal neutral 
axis. The effects of lateral bending, or bending about the 
vertical neutral axis, is neglected. However, for ships 
travelling oblique to the waves, the lateral loads reach a 
magnitude equal to or greater than the vertical bending 
moment. The combined effect of lateral and longitudinal 
bending may be critical for some ship types. 

A requisite for an integrated probabilistic design 
procedure is a more precise knowledge of the loads acting on 



*The numbers in brackets refer co references listed 
in the Appendix 
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a ship. This paper represents a small portion of the total 
effort to completely define the structural loads experienced 
by a ship in a seaway. It is an investigation of the magnitude 
of combined horizontal and vertical bending in irregular seas. 
Due to their relative phase difference, the amplitudes of the 
vertical and horizontal moments generally are not acting at 
the same time. 

In the analysis, the loads are treated as random 
variables. Consideration of the total longitudinal stress at 
the deck leads to the definition of an equivalent combined 
bending moment. Two methods are developed to generate the 
statistical parameters of the combined load so that the usual 
statistical methods to predict long-term extreme values of the 
load can be applied. An example application is then made for 
a large tanker which had been previously instrumented so that 
full scale data for the vertical bending stresses experienced 
in service are available for purposes of comparison. Such a 
comparison serves to check yet another time the validity of 
computational techniques to simulate ship motion and loads in 
waves, and in particular the MIT 5-D seakeeping program. The 
predicted moment responses are compared and analyzed for trends. 



Chapter II 



PROCEDURE 

II. 1 Determination of the Equivalent Combined Wave Bending 
Moment 

A ship travelling at an angle to the waves will be 
subjected to unsymmetrical bending such that the stress at 
the deck is not uniform but assumes a maximum value at one 
edge, as shown in figure II. 1. This primary stress may be 
thought of as the combined effect of vertical and horizontal 
bending and, to a much lesser extent, torsion. Previous 
studies [15, 16] have established that torsional moments are 
relatively small, never exceeding 10% of the vertical moment 
in any given condition. In this study the effects due to 
torsion are considered negligible, but it is recognized that 
torsion loading may be critical for some ships, as those 
with large hatch openings C17]. The contribution to the 
ship primary stress may be included, if so desired, in a 
manner analogous to the method described here [18]. 

The stress at the deck edge is given by 
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The Longitudinal Stress Distribution 
on a Ship Hull Girder 
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where 



Q 

II 


combined, or total, stress 


Q 
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II 


stress due to vertical bending 




stress due to horizontal bending 



In terms of moments, one can write 



where 


^ ^ 
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vertical bending moment 
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horizontal bending moment 
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— = vertical section modulus 
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■5^= horizontal section modulus. 
Y 


multip! 


Lying both sides of the equation by k^, one gets 
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If one defines an equivalent combined bending moment 

K 

Mr = '’tN,- and the section modulus ratio K = 
gets the relation 



then one 



+ KM^ (1) 

It is emphasized here that one cannot simply add 
the maximum values of the vertical and horizontal moments to 
get the maximum combined effect. In general the moments at 
any point along the length of the ship are not acting in 
phase. Furthermore, it may be of interest to specify the 
complete distribution of the primary stress at the deck. 

Within the limits of the beam theory, this is done by getting 
the stress at one other point across the ship's deck. A 
convenient reference point is at the ship's centerline where 
the stress due to lateral bending is zero for most ships due 
to port-and- starboard symmetry, i.e. , at the centerline. 

In the statistical approach to sea loads, the ex- 
citing force is assiamed to be a stationary, narrow-band and 
gaussian process with zero mean over a short period of time 
(e.g. 30 minutes) . The ship is considered to be a linear 
system so that the response is also stationary, narrow-band 
and gaussian with zero mean over the same period of time. 

Then the amplitudes follow the Rayleigh distribution [19, 20]. 
Over a long period of time, the assumption of a stationary 
excitation process cannot be made. However, it was found 



that the amplitudes follow closely the Weibull distribution 
C8, 9, 10]. It is shown that the Rayleigh is a special case 
of the Weibull distribution given by: 

f^(x) = (s/k)(x/k)® ^ exp-(x/k)® X20 (2) 

Fj^(x) = 1 - exp-(x/k)® X>0 (3) 

wher probability density and distrib- 

ution functions respectively, and s and k are parameters. To 
get the Rayleigh distribution, one sets s=2 and k = /E where 
E is the mean square value of the amplitude of the response. 
Furthermore it v/as found that the best-fit Weibull curve is 
approximately the exponential distribution obtained from 
equations (2) and (3) by setting s = 1 and k = m where m is 
the mean value of the response amplitude. The statistics of 
the ship bending moment response, therefore, can be specified 
if the mean square value, or variance, of the process is 
known . 



The usual method to obtain the variance of the re- 
sponse in random seas, attributed to St. Denis and Pierson, 
consists of specifying the energy density spectra of the 
ocean waves and obtaining the ship response amplitude opera- 
tor (RAO) , or the response per unit amplitude of regular 
wave, by model tests or some computational procedure using 
the strip theory approach. From the theory of linear systems, 
the spectral density of the response Cc<n be found from the 
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relation [21]: 



Sj^(w) = |H(to) 1 ^S(to) 



where S (w) = spectral density of the response 

S\ 



H(co) = ship response amplitude operator (RAO) 
S(w) = spectral density of the seaway. 

The area under the response spectrum is the variance 
of the process and represents the sum of the responses to the 
infinitely many sine waves that are thought to compose the 
actual sea. 

Two methods used to obtain the variance of the wave 
bending moments by computer simulation techniques are ex- 
plained below. 

1) First Method 

Equation (1) can be written showing the moments as 
a continuous random function of time, t: 



M^(t) = M^(t) + KMj^(t) 



mm 




\ 
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Using the assumption that the sea can be represented by in- 
finitely many harmonic waves of small amplitudes, continuous 
frequency and random phase, one can write the vertical and 
horizontal moment responses as: 



M^(t) = V^cos (w^t + e^) , and 

M (t) = H cos (w t + 0 ) 

H n n n 



where 



= amplitude of the vertical moment response to the 
nth wave 

= amplitude of the horizontal moment response to 
nth wave 

0)^ = frequency of the nth wave 

= phase angle of the vertical moment with respect 
to a wave with its crest amidship 

0^ = phase angle of the horizontal moment with respect 
to a wave with its crest amidship 



then 



M^(t) 



CO CO 

I, V cos (w t + e ) + K cos(w t + 0 ) 

n-1 n '• n n^ n=l n ^ n n ■' 



or 



M,j.(t) 



on 




+ KH COS (w t + 0 ) 

n ' n n 'J 



which can be rewritten as 
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M„(t) = I 
n=l 



(v cos e + KH cos 6 ) cos w t 
' n n n n' n 



_ ( 



(v sin e + KH sin 6 ) sin oj t 
'• n n n n^ n 



Let 



X = V cos e + KH cos e 
n n n n n 

Y = V sin e + KH sin 6 
n n n n n 

Then, 

Mm(t) = I X cos (0 t + Y sin o) t 
^ I n n n n 

n=l ^ 

Since M^(t) is also a continuous random function of time 
can also be represented by; 



M,n(t) = y Z cos (w t + 6 ) 

T ' ' S n n n'' 

n=i 



A comparison of equations (4) and (5) shows that 



-1 Y 7 7 2 

e = tan n, and Z = X + Y 
n Y~ n n n 

n 



(4) 

, it 



( 5 ) 
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Then 

Z ^ = V ^ + K^H ^ + 2KH V cos(e - 9 ) (6) 

n n n nn^nn'^ 

Equation (6) gives the amplitude of the combined 
moment due to vertical and horizontal bending, with the 
phase difference in the occurence of the component peaks 
taken into consideration. The MIT 5-D seakeeping program 
calculates the ship motions and the vertical and horizontal 
moment responses to regular waves (RAO) . A slight modifica- 
tion of the program to obtain the equivalent combined moment 
response operator was made by use of equation (6). The RAO's 
were then used in the usual power spectral analysis in the 
frequency domain to obtain the variance of the responses in 
irregular seas. Comparison of predicted and experimental re- 
sults have shown generally good agreement [15, 16, 21, 22, 23] 
and have led to widespread acceptance of the superposition 
method of St. Denis and Pierson. 

2) Second Method 

Considering the wave loads as random variables, it is 
generally accepted that the bending moments M^ and M^ are sta- 
tionary, narrow-band, gaussian processes with zero means for 
short periods of time. Then M^ is also stationary, narrow-band 
and gaussian wirh zero mean, i.e., in statistical notations: 

E M^ 



= E M^ + KE Mjj 



0 



Furthermore 



e[m^] = e[(m^ + KMj,)2 



E^M^j = eI^M^J + K^e[m^ 



+ 2KE 



• M 



H 



but the covariance is defined as: 



E [v^h] P ~ "V ) . ^ _ (^H ~ "'h ) _ } 



where = mean values = 0 



= correlation coefficient of the vertical and 
horizontal bending moments 



therefore. 



E 





+ k2e[m2] 



+ 







E 







For a zero-mean process, the mean square is equal to the 



variance . 
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As in the first method, the value of the section 
modulus ratio, K, must be known for a given ship or esti- 
mated during preliminary design. An attempt was made to 
examine the section modulus ratio of several sizes of tankers 
to determine if there was a trend with the growth in the size 
of these ships. As shown in figure II. 2, the points are 
widely scattered with no apparent trend. Since the value of 
K depends on the construction of the ship (framing system, 
scantlings, etc) , perhaps the scatter of points is indica- 
tive of the present uncertainty about both sea loads and the 
strength of ships, and the rather arbitrary way a "safety 
factor" is arrived at. 

Likewise, the value of the correlation coefficient, 
of the vertical and horizontal moments must be deter- 
mined. Implicit in the use of equation (7) is the assumption 
that an average value of already takes into account the 

phase difference of the peaks of and in all conditions 
at sea. 



With this method, the bending moment response ampli- 
tude operators may be determined by experiment or computer 
techniques. Spectral analysis could then be performed to 
obtain the variances of and and finally the combined 

T. 



moment M, 



ection Modulus 




O > I 22 
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II. 2 Long-Term Wave Load Analysis 

After the statistics of the short-term response are 
determined for the range of sea states that make up the 
actual seaway, the expected long-term loads may be predicted. 
Several methods to obtain the long-term distribution of wave 
loads and the prediction of extreme values have been proposed 
C6, 9, 11, 18, 23]. Mansour Cl8] has described procedures 
for the analysis of failure under extreme loads under two 
conditions : 



1. Short-term analysis - when the route of the ship 
is not known and a "severest condition" criterion 
is imposed. 

2. Long-term analysis - If the ship's route is known 
and if that route is more or less permanent. 

In this study, a procedure is discussed to determine 
the value of the parameter k = m of equations (2) and (3) 
for a ship travelling oblique to waves in all sea conditions 
described by the significant wave height and period (h,T) as 
in [24]. The following assumptions are made: 

1. wave headings are uniformly distributed over all 



directions . 
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2. The speed is essentially determined by the state 
of the sea. The speed-significant wave height relation is 
adopted from [14] as follows: 



Speed, V knots 


significant wave height, h ft. 


16 


less than 12.5 


12 


12.5 to 21.5 


6 


greater than 21.5 



The procedure followed is described below: 

1. Obtain the RMS response to all sea states (i.e., 
range of sea spectra) in all conditions of speed 
and wave headings by some analytical method as 
the "MIT 5-D seakeeping program." 

2. For each speed, get the average of the RMS 
response over all wave directions. 

3. Use Bretschneider ' s two - parameter family of sea 
spectra to get the response to any sea described 
by a particular (h,T) pair. This is most easily 
done by scaling the responses to the sea spectra 
of step (1) . We call this matrix of the response 

R with elements r. .. 

ID 

4. For each wave height in [24], determine the 
probability (i.e., relative frequency of occurence) 




i 





18 



of each observed period. We call the resulting 
matrix of probabilities P with elements p^ ^ . 

5. Determine the average RMS response for each 
significant height by the relation 



m. = y r . . X p. . 
1 j ID 



where the subscripts i and j refer to rows and 
columns of a matrix. The resulting column 
matrix is called M. 



6. Determine the probability of occurrence of each 
significant wave height. We call this column 
matrix Q with elements q^. 

7. The mean or expected value of the RMS of the 
response for a particular area is given by 

mean RMS = 7 m. x q. 

r 1 
1 

Then the desired parameter m = 1.25 x (mean RMS 
of the response) in the case of a short-term 
analysis [18]. For the long-term analysis, a 
time factor is additionally introduced to get the 
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weighted mean of the RMS of the response for all 
areas in the ship's route. 

Steps 4-6 is equivalent to getting the probability 
of occurrence of each (h,T) pair in any one area being de- 
scribed. If we call this matrix P* with elements p! then 

ID 

mean RMS = 7 7 r. . x p’. . 

i j ID ID 

The data in [24] have been computerized and could 
be used in the above calculations. 
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II. 3 Details of the Procedure 

The methods for wave load analysis and prediction, as 
explained in the foregoing sections of this chapter, were 
applied in the analysis of a large tanker. The UNIVERSE 
IRELAND was one of five vessels of different sizes instrumented 
for full-scale measurement of stresses experienced in actual 
service. However, the strain gage arrangement on deck was 
such that the stress due to lateral bending was cancelled out, 
so that only vertical bending stresses were recorded. There- 
fore, the results derived for the response of primary interest, 
the combined moment, could not be verified. The description 
and analysis of the results of the 3^-year project sponsored 
by the American Bureau of Shipping (ABS) is reported in [25] 
and provide the basis for comparison with the analytically 
derived values of vertical bending stress at amidship. 

The general characteristics of the UNIVERSE IRELAND 
are listed below: 



Full load displacement, tons 375,811. 

Length overall, ft 1135.17 

Length between perpendiculars, ft ... 1083. 

Breadth, ft 174.87 

Depth, ft 105. 

Design draft (keel) , ft 81-0-5 

Block coefficient (^1) 0.86 

Section modulus, top, in^-ft 566,794. 

Section modulus, side, in^-ft 617,350. 



The ship's lines and other data pertaining to the 
ship necessary to generate the input for the computer analysis 
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were obtained from the ABS. It was necessary to use an 
approximate lightship weight distribution because of an 
apparent discrepancy in the distribution curve provided. 

Since the still water bending moment, at least, is very 
sensitive to the distribution of the weight, a probable 
error was anticipated in the calculated values. For ex- 
ample, the maximum still water bending moment calculated 
by the computer exceeded the value determined by the ABS 
by about 11.5%. 

The MIT 5-D seakeeping program is fully described 
in [26]. Only minor modifications were necessary to do the 
quasi-static combined bending moment analysis. The re- 
sponse to long-crested irregular seas were used in conjunc- 
tion with an auxiliary program to obtain the responses to 
short-crested seas. This program uses the cosine-squared 
spreading function 

s (/M) = j cos^ (^) 

where /<■= angle between the wind direction and the wave 
components. The principal wind directions were given ex- 
actly the same values as the wave-to-course headings. 

The computer was run for ship's speed of 6, 12 and 
16 knots, nine wave headings from 0° to 180° in increments of 
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22.5*^, and nine sea states described by the significant 

height and period. The computer output consisted of the 

root mean square of the response, the significant height 
1 / 1 / 

(H'*) and the average of the 1/lOth highest (H of the 

response measured peak-to-peak for all combinations of the 

three speeds, nine wave directions and nine sea states. 

Since it was desired to compare the values of the vertical, 

horizontal and combined wave moments, output was obtained 

for all three responses. Furthermore, to be as accurate 

as possible, the three responses were calculated for twenty 

stations specified along the length of the ship, even if 

eventually only the responses amidships were utilized in 

the analysis. 

The sea data were carefully chosen to match the 
actual seaway experienced by the ship during the instru- 
mentation period. This was done by obtaining sea data from 
[24] for the areas lying along the service route of the 
UNIVERSE IRELAND, which made regular runs between the Persian 
Gulf and Western Europe by way of the cape of Good Hope [251. 
Although the seakeeping program can accept and use an actual 
sea spectrum or simulate developing or decaying seas by the 
Bretschneider two-parameter spectrum, the Pierson-Moskowitz 
one-parameter family for fully developed seas was chosen. 

Nine sea states with significant heights ranging from 3.6 feet 
to 43.3 feet (corresponding to observed periods of 1.5 seconds 
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to 18.5 seconds) were used. Spectral ordinates for forty 
frequencies were calculated while twenty five points were 
used to describe the response amplitude operators covering 
the wave frequency range 0.1 < ^ < 9.0 with most of the 
points spaced closely around the expected spectral peak at 
^ = 1.0. X is the wavelength and L is the length of the 
ship. The response to regular waves (RAO's) were not 
printed in order to save on the cost of running the computer. 



Chapter III 



RESULTS 

The following general areas of study were explored: 

1. Comparison of the derived results with those 
in other published works on the subject. 

2. Comparison of the variance of the combined 
moment obtained from the two methods described 
in the previous chapter. 

3. Comparison of theoretical calculations for the 
vertical bending stress with full-scale 
measurements . 

4 . Examination of bending moments trends and wave 
load prediction. 

The results of the above studies are presented in 
the following sections. 
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III.l Comparison of Results with Other Published Works 

The calculated value of the root-mean-square (RMS) 
of the bending moment responses are plotted in figures III.l 
to III. 4 as a function of the wave headings. For contrast 
this was done for the responses to long-crested and short- 
crested fully developed seas with significant wave heights 
of 7.4 ft., 19.2 ft. and 43.3 ft. It was intended to show 
the response to a low, moderate and high sea states. Since 
only the value of the horizontal moment in "realistic" seas 
are of interest, the RMS of the response were plotted only 
for short-crested seas. The graphs were prepared for a 
representative ship's speed. The general behavior of the 
response did not change with speed. 

Faltinsen [22] describes the comparison between 
theory and experiments of tne wave-induced horizontal and 
vertical forces and moments of a series 60 hull form with a 
block coefficient of 0.80 and a Froude member of 0.15 moving 
in regular waves with different headings at an angle to the 
direction of the waves. The predictions made by the new 
strip theory developed by Salvesen, Tuck, and Faltinsen, 
which is used in the MIT 5-D seakeeping program proved very 
satisfactory. The theory is described by Salvesen et al, in 
[27]. For vertical bending moments, the theory predicted 
higher values than the experiments for ^ < 0.7 in quartering 
and following waves and lower values for A/L < 0.6 in bow 
waves. This trend may also be observed by comparing the RMS 
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RMS of the Bending Moment Response at Amidship VS Heading 
in Sho r t — C r e s t ed Waves for Moderate Seas 
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values of in short-crested seas for wave headings of 0*^ 
(following seas) and 180° (head seas) in figures III. 2, 3 
and 4. For horizontal bending moments, appreciable diff- 
erence existed between theory and experiments only for beam 
seas . 



One of the early investigations of the effect of 
horizontal moments in oblique seas is described by Numata 
Cl6]. It was established in this study that lateral moments 
reach values that are substantial and that peak bending 
moment exists in the region of the speed for synchronous 
pitching and heaving motions. Wahab El5] shows that the 
importance of the horizontal bending moment increases with 
increasing ship length. For the sea condition used in his 
investigation, the significant midship horizontal moment for 
a ship of 200 meter-length was 46 per cent of the significant 
vertical moment while for a 400 meter-length ship the per- 
centage increased to 69. The concept of effective wavelength, 
L^, described by Numata was helpful to explain the results 
obtained from the computer analysis. However, only one sea 
spectrum was used in his study and some conclusions were 
reached that are not necessarily true for all conditions at 
sea for any given ship. Inasmuch as the RAO of the UNIVERSE 
IRELAND was not recorded by the computer, the figures in 
[15] proved very useful. The ship studied in that work is 
the same as the model described by Faltinsen and has general 
characteristics similar to the UNIVERSE IRELAND. The response 



of the model to regular waves could be assumed to be sub- 
stantially identical in shape and trends. The calculation 
of amidship forces and bending moments in [15] also made 
use of a single spectrum representing a severe sea. 
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III. 2 Comparison of the Calculated Values of the Combined 
Moment Response ~~ ' 

Table III-l is a summary of the results obtained 
for the RMS value of the combined moment response for the 
vessels' service speeds of 16 and 6 knots. The first set 
of values are those calculated from the RAO developed by the 
computer by use of equation (6) of the first method. The 
second set of values were calculated from equation (7) of 
the second method using the value /yjj = 0.32 obtained from 
[18] as reprinted from the published result of a Det Norske 
Veritas Research Report. The reprinted source is shown as 
figure III. 5. The values of represent the correlation 

coefficient of the most probable largest values of the 
vertical and horizontal moment for one ship life of about 
twenty years in the North Atlantic. They were meant to be 
reference values and not for use in design. The difference 
in M^ values is shown for each sea state. It is immediately 
apparent that the results obtained by both methods are not 
very different. The third set of values was calculated also 
by use of equation (7) but using a value of obtained by 

averaging the responses in short-crested seas by the first 
method for all wind-wave directions and for three represent- 
ative spectra for low, moderate and high seas. The mean 
value of p.j„ obtained in this manner was 0.5 32, significantly 
higher than the ISSC value. The difference in the results for 
M^ is also indicated for comparison. 



Comparison of the RMS Response for the Combined 
Moment at Amidships as Calculated by 
the Two Methods for Two Speeds 
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Comparison of the RMS Midship Combined Moment Response 
as Computed by Two Methods 



35 



Series 60, All Headings Included, Q=lo“® 



c^z0.<a 

Ce^o.7 Fn"0 20 

Cgr o.S F,^= o \5 




lOO 2.00 300 L (lrv^') 400 




>oo aoo 300 




\oo zoo 300 \-(^VTV) 400 



Figure III. 5 

Correlation Coefficient Between Vertical Bending Moment 
and Horizontal Bending Moment 
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III. 3 Comparison of the Calculated Vertical Bending Stress 
and Full-scale measurements ; 

Strain gages were installed on the port and starboard 
sides at amidship and several other stations along the length 
of the UNIVERSE IRELAND to measure stresses on the deck. The 
arrangement was such that only vertical bending stresses were 
measured. Therefore, though the combined moment could not be 
verified, a comparison between full scale data and the com- 
puter-derived values was possible. This comparison by 
weather group as specified by wind measurements (Beaufort 
scale) is shown in figure III. 6. A summary is also pre- 
sented in table III. 2. The conversion from bending moment 
to stress is done by the relation: Bending Moment = Stress x 

Section Modulus. The values compared are the average RMS 
values of the peak-to-peak stress for each wind scale. The 
full scale data were read off figure 11 of [25]. The con- 
version from wind to wave was done by the graphical relation 
shown in figure III. 7 C23]. This conversion is summarized 
in table III. 3. 



The use of nine sea spectra with each significant 
wave height enables the calculation of the standard devia- 
tion of the predicted values of the vertical stress. For 
comparison with the measured data, these standard deviations 
had to be modified to obtain values based on wind measure- 
ments. This was accomplished by the relation [23]: 



37 



S 



1 



2 




_2 
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where = standard deviation of the response with respect 

to the wind 

S 2 = standard deviation of the response with respect 
to waves 

= standard deviation of the wave height with re- 
spect to the wind, as shown in figure III. 8 
from [23 J. 

Ah = wave height increment 

0 = average slope of the response vs. significant 
height plotted in figure III. 9. 

For comparison, table III. 4 of the standard deviations for 
all the cases studied is drawn up. 
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Comparison of Predicted Average RMS Peak-to-Trough 
Vertical Bending Stress With Measured Values 
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Relationship Between Significant 
and Wind Speed from Various 



Wave Height 
Sources 
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Figure III. 8 

Wind- Wave Height Standard Deviations 
from Various Sources 



Summary of Wind-to-Wave Conversion 
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III. 4 Examination of Bending Moment Trends and Wave Load 
Prediction 

Since the traditional primary strength calculations use 
the vertical bending moment in a long-crested sea of arbirtary 
height as a basis for ship load, it is of interest to compare 
the values of the vertical bending response in long-crested head 
seas with the combined moment in short-crested seas for the 
higher range of the wave significant heights. Figures III. 10 
and III. 11 are drawn for sea spectra representing fully developed 
moderate and severe seas, respectively moderate seas, the com- 
bined moment for head seas exceeded the long-crested vertical 
moment by about thirty sever per cent (37%) and in severe seas, 
the difference was sixteen per cent (16%) . 

Similarly, it was deemed desirable to compare the rel- 
ative values of the vertical and horizontal moments for all 
conditions at sea. This is done in figures III 2, 4 and 6 for 
the responses in short-crested seas. One can see that the 
horizontal moment is not small compared to vertical moment; 
and in low to moderate seas, the former is at least as great as 
the latter. The above mentioned comparison of the response in 
long-crested and short-crested head seas reflect the relative 
contributions of the vertical and horizontal moments to the 
total load. 

Table III. 5 illustrates the change with speed in the 
bending moment response for low, moderate and high seas. The 
tabulated values are the mean RMS peak-to-mean responses in 
short-crested seas averaged over all wind-wave directions. 
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RMS Single Amplitude 
Combined Bending Moment 
Hp(f t-tons) 


163295.8 

159227.2 

159172.3 


829528.0 

828541.7 

820420.8 


1651275.8 
1657771.0 

1662724.8 


RMS Single Amplitude 
Horizontal Bending Moment 
M^(f t-tons) 


1 

i 112629.9 

1 

110631.4 

110136.3 


490273.9 

506827.1 

510401.6 


811011.7 

845107.0 

870121.2 


RMS Single Amplitude 
Vertical Bending Moment 
M^(f t-tons) 


76668.6 
72738.8 

70636.6 


473399.7 

463361.7 
453427.6 


1307839.4 

1146911.1 

1133966.6 


speed 

(knots) 


16 

12 

6 


16 

12 

6 


16 

12 

6 


Significant 
wave height 
(ft.) 


II 


h = 19.2 


h = 43.3 
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In order to illustrate the procedure to obtain the 
parameter k for the long-term distribution of the wave re- 
sponse, a "short-term" analysis was done for the worst 
condition in the ship's route (area 7 in [24]). This is the 
North Atlantic in the months of December to February and the 
tabulated data of [24] were used. The results obtained are 
given as tables III. 6 to III. 9 for the combined moment. For 
this analysis, the result is 

mean RMS = 367820 ft-tons 

Since this is a "short-term" analysis, one can then say 
m = 1.25 (mean RMS) = 459775 ft. tons. 

The equivalent values for the vertical and horizon- 
tal moments were also computed. These are: 

mean RMS response, M^ = 197307 ft-tons 

mean RMS response, M„ = 231043 ft-tons 

Note that the expected value for the horizontal moment is 
larger because the ship encounters more low-to-moderate 
seas in any length of time than severe seas. 



The mean value of the RMS of the response was also 
calculated for all the areas in the route of the ship, i.e. 
a "long-term" analysis was conducted. Both and were 
calculated and the results are summarized in table III. 10. 
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Table III. 6 

Mean RMS of the Combined Moment Response in Short-Crested Seas 
Averaged Over Uniformly Distributed Wind-Wave Directions 



Velocity 

(knots) 


Significant 
Wave Height 
(ft) 


RMS My 
(f t-tons) 


RMS Mjj 
(ft- tons) 


RMS 

(ft- tons) 




3.6 


13386.3 


21411.2 


26762.9 




7.4 


54212.8 


79641.4 


115467.5 




10.7 


115913.9 


159110.8 


240620.1 




14.6 


212240.4 


253534.2 


405516.0 


V=16 


19.2 


334744.1 


346676.0 


586564.9 




24.3 


463645.8 


426305.6 


756493.1 




30.0 


591299.8 


488621.4 


911865.7 




36.4 


712415.3 


537815.8 


1050905.6 




43.3 


924782.1 


573471.8 


1,167,628.3 




3.6 


11689.7 


19464.6 


24197.6 




7.4 


51434.1 


78228.2 


112,590.6 




10.7 


112157.7 


160,233.2 


238895.4 




14 . 6 


206043.0 


258953.5 


403795.1 


V=12 


19.2 


327646.2 


358380.9 


585867.4 




24.3 


456173.6 


442364.9 


757601.8 




30.0 


583590.6 


509395.5 


914721.5 




36.4 


705103.0 


561731.7 


1055359.7 




43.3 


810988.6 


597580.8 


1172221.1 




3.6 


9831.6 


17268.6 


22238.0 




7.4 


49947.6 


77878.1 


112551.8 




10.7 


110104.1 


157927.6 


236043.6 




14.6 


201175.2 


258401.3 


398461.4 


V=6 


19.2 


320621.7 


360908.4 


580125.1 




24.3 


449758.8 


449322.7 


756008.3 




30.0 


696331.0 


518255.7 


913820.7 




36.4 


696331.0 


573820.7 


1056466 . 6 




43.3 


801856.7 


625268.6 


1175724.0 



The Response Matrix R 
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Table 7 (continued) 
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Table III. 8 (continued) 
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Short-Term** Analysis of Area 7 for Months of December-February . 
Weighted RMS Values for the Combined Moment Response 
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Table III -9 (continued) 
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Table III. 10 

Mean RMS Value of the Response for all Areas in the 
Route of the UNIVERSE IRELAND 



Area 


E[RMS]of M^ Response 


E[RMS]of Response 


Time Factor 


7 


154254 


299294 


2 


11 


126213 


246724 


2 


18 


109398 


218546 


1 


20 


66961 


138544 


1 


23 


110971 


220161 


1 


28 


80454 


170798 


1 


29 


92992 


189706 


1 


34 


108754 


203516 


1 


41 


146576 


282585 


2 


42 


123456 


246480 


2 


45 


161230 


260578 


2 



Time Weighted Average : 
(RMS)^ = 124,562 ft-tons 

(RMS)^ = 517,450 ft-tons 



“v 

y 



= 155,702 ft-tons 
= 646,813 ft-tons 

= 4.15 



Chapter IV 



DISCUSSION OF RESULTS 

IV. 1 Evaluation of the analytical Results 

The results obtained in this study agree with the 
previous findings in other studies that the lateral moment 
can be very significant for ships in waves. Some definite 
trends of the moment responses can be observed from figures 
III.l to III. 4. From the results in short-crested seas 
shown in figures III. 2, 3 and 4, one sees that horizontal 
and vertical moments change very little with wave direction 
in moderate seas. On the other hand, there is a broad peak 
away from the head and following seas for all moment responses 
in low sea states. In severe seas, the maxim\m response for 
and are in the head and following seas while stays 
almost constant regardless of wave heading. For the long 
crested seas of figure III.l, one can see that the vertical 
moment peaks at about 60° for low to moderate sea states and 
assumes extreme values at head and following seas in high 
seas. The equivalent combined moment has two peaks in all 
cases. The peaks occur in wave headings near 60° and 120° 
with a shift away from the beam seas case as the sea grows 
more severe . 
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Vcriical wave bending moment 
amplitude in regular bow waves 



Figure IV . la 




Vertical wave bending moment 
amplitude in regular beam and 
quartering waves 

Figure IV . lb 





Horizontal wave bending moment 
amplitude in legular bow waves 



Horizontal wave bending moment 
amplitude in regular beam and 
quartering w'aves 



Figure IV. 2a 



Figure IV. 2b 
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One observes from figure III. 2 and III. 3 that the 
horizontal moments exceed the vertical moments in almost all 
headings for low and moderate seas. However, in severe seas, 
for which structural design is based, one sees the customary 
relation of the vertical moment exceeding the horizontal 
moment, especially in head and following seas. The reason 
for this phenomenon is that lateral bending is more respon- 
sive to higher frequency excitation than vertical bending. 

One can see from figures IV. 1 and IV. 2 that the spectral peak 
of the lateral moment response for all wave headings are gen- 
erally situated at higher frequencies than the corresponding 
response curves for the vertical moments. One can note fur- 
ther that the spectral ordinates of the RAO for wave headings 
70°, 110° and 130° are higher for the lateral moment than the 
vertical moment. Therefore, when the energy of the exciting 
waves are concentrated in the higher frequency range (i.e., 
in low sea states) , the lateral moment exceed the vertical 
moment. In high sea states, the reverse trend is true. 

For a similar reason the peaks of the response curves 
for the vertical moment occur away from head and following 
seas in low sea states. The shift in the location of the 
spectral peak of the RAO's in the direction of higher cir- 
cular frequency, u, can be explained by considering the 
effective wavelength in oblique seas. It is generally known 
that the wavelength nearly equal to the ship length causes 
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V 




Figure IV . 3 

The Effective Wavelength in Oblique Waves 
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the maximum midship bending moments. In oblique seas the 
ship senses the effective wavelength, L^, as shown in 
figure IV. 3. The relation with the actual wavelength is 

= L/cos ^ 

where = effective wavelength 

L = actual wavelength 
^ = heading angle 

As ^ approaches 90° from either 0° or 180°, the actual 
wavelength decreases while the value of stays fixed. 
Therefore the peak frequency of the response goes up, thereby 
shifting the RAO spectrum nearer the sea spectral energy 
density peak in low sea states. 

The theoretical calculations predict greater vertical 
moment response in following seas than in head seas, contrary 
to experimental results. Wahab's results [15], shown as 
figure IV. 4, confirm this apparent discrepancy betv/een the 
theory and experiments. 

Finally, one notes that the combined moment in 
"realistic" seas is significantly higher than either of the 
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components in all cases but, because of the relative phase 
angles, is not equal to the sum of their individual peaks. 
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IV. 2 Comparison of the Values of the Combined Moment 

Table III.l shows the comparison of the response 
RMS value of the combined moment, as obtained by the two 

methods described in chapter two. It can be readily seen 
that the results are not significantly different, the aver- 
age for all sea spectra being on the order of eight per cent 
(8%) . This confirms the statement that the average value of 
the correlation coefficient, already considers the phase 

difference of the vertical and horizontal bending moment 
peaks. As would be expected, the third set of values for 
were closer to the computer values. It is realized that the 
two values of that were used were derived in a different 

manner and perhaps should not be directly compared. In de- 
riving the second value, it was found that the coefficient 
varies strongly with the state of the sea and the wave di- 
rection and therefore that no single value will be "true" for 
all conditions at sea. If this method is to be adopted, 
there is need to study the correlation coefficient further. 
There is a dearth of published material on the subject. 



The calculated value of the combined moment also de- 
pends on the relative phase angle of the vertical and hori- 
zontal moments. Therefore, it is necessary that the pre- 
dicted values be at least as accurate as the predictions for 
the vertical and horizontal moment responses. Comparison of 
theoretical and experimental phase angles shows that the 
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agreement is in general quite good [22]. The maximum dis- 
crepancies occurred for wave headings of 50° and 110° for 
the wave lengths X/L = 1.2 and ^ = 1.5. (Figures IV. 6 and IV. 7) 

In view of the above mentioned considerations it 
would seem to be indicated that both method 1 and method 2 
to obtain the combined moment are equally "good". The first 
method is probably more convenient for use in analytical 
studies . 
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Figure IV . 6 



Phase Angle Between Horizontal and Vertical 
Wave Bending Moment at Midship 
for Series 60, Block 0.80 
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Figure IV . 7 



Phase Angle Between Horizontal and Vertical 
Wave Bending Moment at Midship 
for Series 60, Block 0.80 
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IV. 3 Comparison Study of the Vertical Bending Stress 
Data on the full-scale measurement of vertical 
bending stresses on the deck of the UNIVERSE IRELAND can be 
found in [253. A comparison between calculated and measured 
values was done by plotting the average RMS peak-to-peak 
amidship bending stresses. The analytical values were ob- 
tained by calculating the response to each of nine sea 
spectra corresponding to the observed periods in the tabu- 
lated data of [243. The significant wave heights are the 
wind-to-wave conversion shown in table III. 3. The average 
RMS peak-to-peak response is the weighted mean of the cal- 
culated response to each sea specified by the significant 
height. This is best described by the relation 

V h 

mean RMS = 2. ’ P 

i=l ^ 



where 

h 

Ey = RMS peak-to-peak bending stress 

p(T) = relative frequency of occurrence of each 
observed period for a given significant 
wave height. 
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The wind-wave relation as described by the Pierson- 
Moskowitz random and ITTC 1966 line, the corrected Roll curve 
for the North Atlantic and the corrected Yamanouchi curve for 
the North Pacific were examined. The conversion is summar- 
ized in table III. 3. 

The use of the Pierson-Moskowitz Random line pro- 
duced the worst case in the comparison study. While the 
v/ind-to-wave relation is intended to convert wind measure- 
ments to significant wave heights for seas in random stages 
of partial development, there is a bias in favor of seas that 
are near full development [25]. Consequently, the values 
based on this wind-to-wave conversion, coupled vzith the known 
tendency of the theory to over-predict [21], are injch higher 
than the measured stresses. The Roll and Yamanouchi curves 
were based on observed wave heights and corrected for signif- 
icant heights. They are believed to be better representation 
of an actual sea. Since the ship route is neither entirely 
in the North Atlantic nor in the North Pacific, one sees fairly 
good agreement with results based on both the Roll and 
Yamanouchi relations for the lower half of the sea states as 
described by the beaufort scale. It can be seen that theoret- 
ical results over-predicted the average values of the RMS 
stress, specially in the higher sea states. There are several 
probable causes of this: 
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1. Since the largest divergence of compared values are 
consistently in the more severe seas, one is led 

to believe that the assumption of linear motions 
for the ship, in particular resonance in roll, 
exaggerates the analytical values. Obviously, 
heave and pitch are large in heavy seas. The 
divergence caused by roll resonance gets bigger 
as the energy of the waves are concentrated more 
and more in the lower frequencies. Faltinsen [22] 
used a first approximation of roll at resonance 
as input values. It is also to be noted that 
the ship was not provided with bilge keels in 
the computer analysis. It is not known if the 
acutal ship has these dampeners . 

2. Loukakis [21] suggested that deflections of the 
hull can relieve some bending stresses, so that 
even if bending moment predictions are accurate 
the resulting conversion into stresses may result 
in over-prediction. 

3. The inexact weight load distribution used in the 
analysis could have caused some of the over predic- 
tion. A qualitative measure of this effect can be 
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gained by the excess in value in the calmer seas 
when the contribution of inertia effects to the 
dynamic loads are relatively greater. 

4 . Exaggeration of the reported occurrence of the 
higher period waves. In high seas, when the ship 
is pitching and heaving, it is understandable 
that such errors could be made by untrained ob- 
servers. In this study, this effect would be 
very minor . 

5. Inaccuracy in the measured data due to insuffi- 
cient number of records in extreme v;eather. 

For direct comparison of the scatter of the derived 
values of the stress with the measured values, it v;as necess- 
ary to determine the standard deviations with respect to 
wind measurements. One can see from table III. 4 that the 
calculated response to the nine sea spectra with the same 
significant wave height were more widely scattered than the 
measured RMS peak-to-peak values. 
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IV, 4 Bending Moment Trends and Wave Load Prediction 

One can observe from figures III. 10 and III. 11 that 
the combined moment in short-crested seas are considerably 
higher than the vertical moment in head seas. The difference 
in severe seas is in the order of 15-20 per cent. Therefore, 
at least for this ship, calculations based on vertical bend- 
ing only could lead to disastrous results. 

Comparison of the horizontal and vertical bending 
moments in short-crested seas have been partially discussed 
in section IV. 1. It was explained that in low-to-moderate 
seas the lateral bending response to wave excitation is 
greater than vertical bending. One consequence of this is 
its possible significance in fatigue failures, .inother is 
that over any length of time the expected wave bending load 
will be greater for lateral moments than for vertical moments. 
This is so because the ship will encounter more low-to-moderate 
seas than high seas in its lifetime. The result for the 
short-term analysis for the North Atlantic confirms this 
statement. The statistics reported in [24] were collected 
for a period of seven years. 

Table III. 5 shows that in short-crested seas the 
bending moment varies only slightly with speed. This state- 
ment holds true for all the moment responses, with 
having the greatest variation for all three sea states shown. 



This indicates that reducing speed in a seaway reduces 
dynamic stresses due to slamming, etc. but will not 
significantly reduce the primary stresses due to wave 
bending loads. 



Chapter V 



The 


CONCLUSIONS AND RECOMMENDATIONS 
subject of the combined effect of horizontal and 



vertical wave moments has been investigated. A number of 



conclusions 


can be drawn from the results of the study. 


1. 


The results obtained by the MIT Seakeeping 
Program agreed well with published results of 
theoretical and experimental investigations 
of ship motions and loads in a seaway. The 
following bending moment trends were reserved: 


i. 


The horizontal moments experienced by the UNIVERSE 
IRELAND are not negligible in comparison with the 
vertical moments in almost all conditions at sea. 
In general, the relative significance of the lat- 
eral moments increase with the size of the ships 
Cl5]. In low sea states, the lateral moment 
exceed vertical moments . 


ii . 


The variation with wave headings of the moment 
response seem to depend on the state of the sea. 
For long-crested seas, the maximum vertical 
moment occurs at a wave heading near 60° for low 
to moderate seas and at 0° and 180° in severe 
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seas. The horizontal and combined moments, on 
the other hand, exhibit two peaks near 60° and 
120° in all sea conditions. For short-crested 
waves, all moment responses show a broad peak 
away from 0° and 180° in low seas. The vertical 
and horizontal moments flatten out in moderate 
seas, and the vertical and combined moment 
assume their biggest values at 0° and 180° in 
severe seas. 

The above behavior of the moments with respect to 
wave headings can be explained by the character of 
the RAO. 

2. The values of the combined moment obtained by the 
two methods described in this paper were not very 
different. Therefore, both methods can be con- 
sidered valid. The first method involves develop- 
ing the RAO for the equivalent combined moment by 
a strip-theory based computational method and per- 
forming the usual power spectral analysis. Invest- 
igations on the correlation between analytical and 
experimental values have led to the general accep- 
tance of the method first introduced by St. Denis 
and Pierson. However, it is felt that further 
investigations and refinements are necessary be- 
fore the analytically derived values can be used 
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in actual design. The second method makes use 
of the correlation coefficient, to take into 

account the phase difference of the vertical and 
horizontal moments. There is also further need 
to investigate this important statistical 
parameter. 

3. The comparison study of the predicted and measured 
vertical bending stresses showed good agreement 

in the lower half of the weather spectra en- 
countered in service. In higher seas, the excess- 
ive predicted values are believed to be caused by 
non-linear motions, contrary to the assumption of 
the theory. No doubt the addition of bilge keels 
in the analysis of the ship would have reduced the 
errors . 

For the ship's route, both the Roll and Yamanouchi 
corrected wind-to-wave relations seemed to be valid 
representations of the seaway. In view of the ex- 
pected higher- than-actual values by theory, Roll's 
curve is deemed the better approximation of the 
service route if used in conjuction with the tabu- 
lated data of [24], as was done in this study. 

4. It was found that there was a significant diff- 
erence between the values of the vertical moment 
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response in long-crested head seas and the 
combined moment response in short-crested seas. 
Again, it is stated that this may not be true 
for other ships, particularly the smaller ones. 

For large tankers with the general character- 
istics of the UNIVERSE IRELAND, the contribution 
of the horizontal moment to the longitudinal 
stresses must be considered in the calculation 
of the primary strength. 

Speed is not a significant factor in the magnitude 
of the wave loads acting on ships. 

It remains to be said that the absolute values of the 
combined moment, or the total longitudinal stress, as pre- 
dicted by the methods in this study could not be compared 
with experimental results. However, the moment trends were 
studied and found to be reasonable and, on the whole, in 
conformity with other researchers' results. Therefore, one 
has confidence on these derived values. The theoretical vert- 
ical bending stress predictions could have been closer to 
measurements but, again, one has firm ideas about the pro- 
bable reasons for the divergence and therefore one is not 
completely dissatisfied. 

The procedure followed in this study is fairly in- 
volved and the cost is not cheap. For further analysis 



of ships in related studies, it is recommended that the RAO 
or an approximation, be obtained in the cheapest manner. 
Hopefully, with some anticipation of the expected results, 
more realistic values can be obtained. 
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Appendix A 



SAMPLE COMPUTER OUTPUT DATA 



I. Irregular Wave Results - Long Crested Unidirectional Seas 





Ship Speed = 


12 kts. Heading 


Angle = 22.5 


degrees 




Significant Wave Height = 19.2 


feet 






Peak Spectral 


Frequency = 0.5197 




Motions 


; at Origin: 










RMS 


H(l/3) 


H(l/10) 


H(l/1000) 


Heave 


0.5974 


2.3897 


3.0469 


4.6003 


Pitch 


. 2456 


0.9822 


1.2524 


1.8908 


Sway 


. 2137 


0.8549 


1.0901 


1.6458 


Roll 


.4386 


1.7546 


2.2371 


3.3775 


Yaw 


.1301 


0.5204 


0.6635 


1.0017 



Dynamic 

Sta 


Loadings : 

Vertical 

RMS 


B.M. 

H(l/3) 


H(l/10) 


H(l/1000) 


1 


0.0 


0.0 


0.0 


0.0 


2 


51.1 


204.4 


260.7 


393.6 


3 


1454.4 


5817.9 


7417.9 


11199.6 


4 


8894.8 


35579.2 


45363.5 


68490.0 


5 


41643.6 


166574.5 


212382.6 


320656.0 


6 


89928.7 


359715.0 


458636.6 


692451.3 


7 


144749.5 


578998.0 


738222.5 


1114571.0 


8 


199801.3 


799205.5 


1018987.0 


1538470.0 


9 


251274.3 


1005097.2 


1281499.0 


1934812.0 


10 


296087.0 


1184348.0 


1510043.0 


2279869.0 


11 


329147-6 


1316590.0 


1678653.0 


2534437.0 


12 


347970.3 


1391881.0 


1774648.0 


2679371.0 


13 


353323.6 


1413294.0 


1801950.0 


2720591.0 


14 


344392.0 


1377568.0 


1756399.0 


2651818.0 


15 


321096.7 


1284387.0 


1637593.0 


2472444.0 


16 


284297.1 


1137188.0 


1449915.0 


2189088.0 


17 


180498.0 


721992.2 


920540.1 


1389835.0 


18 


65363.6 


261454.4 


333354.4 


503299 . 8 
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M 



Vertical B.M. (cent.) 



Sta. 


RMS 


H(l/3) 


H(l/10) 


H(l/1000) 


19 


24263.0 


97052.3 


123741.6 


186825.7 


20 


3266.2 


13065.0 


16657.9 


25150.2 


21 


3990.7 


15963.1 


20353.0 


30729.0 


22 


3683.9 


14735.9 


18788.3 


28366.7 



Sta. 


Horizontal 

RMS 


B.M. 

H(l/3) 


H(l/10) 


H(l/1000) 


1 


0.0 


0.0 


0.0 


0.0 


2 


4568.3 


18273.5 


23298.8 


35176.6 


3 


15469.4 


61877.6 


78894.0 


119114.4 


4 


24915.5 


99662.0 


127069.1 


191849.4 


5 


52746.4 


210985.7 


269006.8 


406147.5 


6 


89010.3 


356041.5 


453952.9 


685379.8 


7 


129957.7 


519831.0 


662784.5 


1000674.6 


8 


170987.7 


683951.0 


872037.5 


1316605.0 


9 


208736.5 


834946.2 


1064556.0 


' ~ “271.0 


10 


238483.0 


953932.0 


1216263.0 


lo^.^l9.0 


11 


258256.3 


1033025.5 


1317107.0 


• 1988574.0 


12 


268522.5 


1074090.0 


1369464.0 


2067623.0 


13 


269648.6 


1078594.0 


1375208.0 


2076294.0 


14 


260647.3 


1042589.5 


1329301.0 


2006984.0 


15 


242403.2 


969613.0 


1236256.0 


1866504.0 


16 


214430.1 


857720.5 


1093593.0 


1651111.0 


17 


132947.2 


531789.0 


678031.0 


1023693.7 


18 


48485.1 


193940.3 


247274.0 


373335.2 


19 


22369.9 


89479.8 


114086.7 


172248.6 


20 


6595.0 


26380.2 


33634.7 


50781.8 


21 


456.7 


1826.8 


2329.2 


3516.7 


22 


118.8 


475.5 


606.2 


915.3 





Combined B.M. 








Sta. 


RI4S 


H(l/3) 


H(l/]0) 


H(l/1000) 


1 


0.0 


0.0 


0.0 


0.0 


2 


4177.5 


16710.1 


21305.4 


32167.0 


3 


15031.4 


60125.7 


76660.3 


115742.1 


4 


30139.6 


120558.4 


153712.0 


232075.0 


5 


86181.8 


344727.5 


439527.5 


663600.3 



Combined B.M. (cont.) 



Sta. 


RMS 


H(l/3) 


H(l/10) 


H(l/1000) 


6 


164735.4 


658941.7 


840150.7 


1268462.0 


7 


253106.8 


1012427.2 


1290844.0 


1948922.0 


8 


342644.4 


1370577.0 


1747486.0 


2638362.0 


9 


424373.1 


1697492.0 


2164303.0 


3267672.0 


10 


492267.6 


1969070.0 


2510565.0 


3790460.0 


11 


540433.3 


2161733.0 


2756210.0 


4161336.0 


12 


566890.0 


2267560.0 


2891139.0 


4365052.0 


13 


571413.7 


2285655.0 


2914210.0 


4399885.0 


14 


554399.9 


2213199.0 


2821829.0 


4260409.0 


15 


513669.5 


2054678.0 


2619714.0 


3955255.0 


16 


452582.7 


1810331.0 


2308172.0 


3484887.0 


17 


281527.8 


1126111.0 


1435792.0 


2167764.0 


18 


99230.1 


396920.7 


506073.9 


764072.3 


19 


37856.3 


151425.4 


193067.4 


291494.0 


20 


6591.3 


27805.4 


35451.8 


53525.4 


21 


4085.7 


16342.9 


20837.2 


31460.1 


22 


3671.7 


14687.0 


18725.9 


28272.5 
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II, Irregular Wave Results - Short-crested Multi-Directional 
Seas 



Ship Speed = 12 kts. 

Significant Wave Height = 43.3 feet 
Peak Spectral Frequency = 0.3461 
Principal Wind Direction = 22.5 degrees 
Motions at Origin: 





RMS 


H(l/3) 


H(l/10) 


H(l/1000) 


Heave 


6.2544 


25.0175 


31. 8974 


48.1588 


Pitch 


1.1184 


4.4734 


5.7036 


8.6113 


Sway 


3.5923 


14.3690 


18.3205 


27.6603 


Roll 


4.1180 


16.4720 


21.0018 


31.7086 


Yaw 


0.6277 


2.5109 


3.2014 


4.8334 



Dynamic 

Sta 


Loadings : 
Vertical B. 
RMS 


M. 

H(l/3) 


H(l/10) 


H( 1/10 00) 


9 


660410.1 


2641640.0 


3368092.0 


5085158.0 


10 


770714.8 


3082859.0 


3930646.0 


5934504.0 


11 


852958.4 


3411833.0 


4350088.0 


6567779.0 


12 


902672.4 


3610689.0 


4603629.0 


6950577.0 


13 


916354.1 


3665416.0 


4673406.0 


7055927.0 


14 


891059.9 


3564239.0 


4544406.0 


6861161.0 


15 


827132.9 


3308531.0 


4219378.0 


6368923.0 


16 


728960.0 


2915840.0 


3717696.0 


5612991.0 


17 


456989.5 


1827958.0 


2330646.0 


3518819.0 



Sta 


Horizontal 

RMS 


B.M. 

H(l/3) 


H(l/10) 


H(l/1000) 


9 


440864.0 


1763456.0 


2248406.0 


3394653.0 


10 


518902.5 


2075610.0 


2646403.0 


3995549.0 


11 


575644.3 


2302577.0 


2935786.0 


4432461.0 


12 


606.50.0 


2424600.0 


3091365.0 


4667354.0 


13 


609593.2 


2438373.0 


3108925.0 


4693867.0 


14 


584602.6 


2338410.0 


2981473.0 


4501440.0 


15 


535193.7 


2140775.0 


2729488..0 


4120991.0 


16 


463863.0 


1855452.0 


2365701.0 


3571745.0 


17 


273484.4 


1093937.0 


1394770.0 


2105830.0 



Combined B.M 



Sta 


RMS 


H(l/3) 


H(l/10) 


H(l/1000) 


9 


958354.0 


3833416.0 


4887606.0 


7379326.0 


10 


1121474.0 


4485896.0 


5719517.0 


8635349.0 


11 


1242414.0 


4969656.0 


6336311.0 


9566587.0 


12 


1310208.0 


5240832.0 


6682061.0 


10088601.0 


13 


1320890.0 


5283560.0 


6736539.0 


10170852.0 


14 


1273295.0 


5093180.0 


6493804.0 


9804371.0 


15 


1170845.0 


4683380.0 


5971309.0 


9015506.0 


16 


1020343.5 


4081374.0 


5203752.0 


7856645.0 


17 


616199.8 


2464799.0 


3142619.0 


4744738.0 



Appendix B 



CALCULATIONS FOR AN AVERAGE VALUE OF 

For comparison with the value of the correlation 
coefficient of = 0.32 as obtained from figure III. 5, 

an average of was calculated from the moment responses 

computed by equation (6) of the first method. 

^ e[m^, 2] - E[n^2j . k2e[mjj2 

VH 

2K E[M^2j 



As the long-term extreme wave moments were not ex- 
trapolated, the mean value was obtained by averaging 

over wind/wave directions and for representative sea states . 
The results of the calculation are summarized below: 
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V = 16 knots 

Short-crested seas responses 



wind/wave direction 


h=7.4 ft. 

Ah 


h=24.3 ft. 

Ah 


h=43.3 ft. 

Ah 


0° 


0.771 


.692 


.590 


22.5° 


.799 


.700 


.598 


45.° 


.633 


. 697 


.603 


67.5° 


.743 


.650 


^ 

:> 


90° 0 


. 676 


.555 


. / / 4 


112.5 


.587 


.441 


.367 


1350 


.460 


.384 


.334 


157.5° 


. 293 


.414 


. 368 


180° 


.214 


.468 


.395 


average 


0.575 


0.556 


0.477 



mean for low, moderate and high seas 
/^VH = 0.536 



Appendix C 



CALCULATION OF THE MEAN RMS PEAK-TO-PEAK 
VERTICAL BENDING STRESS AND THE 
STANDARD DEVIATION WITH RESPECT TO WIND 

The calculated bending moment was converted into stress 
by the relation: 

Stress = Moment -i- Section Modulus 

where the section modulus of the ship at the deck is given as 
566,794 in^ - ft [25]. 

The computer output gave the RMS of the response. To 
get the peak-to-peak RMS value, the data was multiplied by 
the factor 2/2. 

The mean value of the response to a sea described by 
its significant height is given by: 

9 

^ = I p(T^) 
i=l 

where m^ = mean RMS of the peak-to-peak response to a sea 
specified by the significant height. 

V. = RMS of the peak-to-peak response to the ith 
^ sea spectra 
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= probability of occurrence of the ith sea spectra. 

is approximated by the relative frequency of 
occurrence of the different observed wave periods for a 
given wave height in the tabulated data of [24]. 

Average values for all the areas in the ship's route 
were obtained. The calculation is illustrated in table A.l 
for the one case of h = 13.2 feet (wave height code 05 in 
[24]) . 



The standard deviation of the vertical bending stress 
with respect to wave measurements is calculated from the 
relation 



S 



2 

2 





p(Ti) 



where S 2 = standard deviation of the bending stress with 
respect to waves. 



For direct comparison with experiments based on 
wind measurements, the following relation is used: 



s 



2 



1 




12 



2 

tan e 



[23] 



Sample Calculations for Table A.l Wave Height Code = 05 (h=13.2 ft) 
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<r 

o 

VO 



o 

m 



VO 

vO 



CO 

CN 
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m 
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CN 
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VO 



o 
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CN 

m 

o 
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m 

VO 
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00 

in 

o 



vO 

CN 

00 



CO 



CO 



00 

o 

o 

CO 



o 

CO 

CN 
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o^ 

VO 

m 

m 

CNJ 



o 

00 

vO 

m 

CN 



Ov 

00 

<r 

CM 



in 



CM 

CM 



VO 

00 

CM 

VO 

00 



CO 

VO 

CO 

00 

CO 



vO 

o^ 

CO 



o 

CO 

CO 



CO 

CM 






CM 

CO 



00 

VO 



VO 

o^ 

rH 

VO 

rH 





00 






CO 








<r 




o 




CM 


m 






CO 




00 


o 


rH 




m 




rH 




OV 










vO 


CM 










rH 


VO 








00 


rH 


00 








VO 


O 


CM 








CO 
















Q) 




CO 






<D 


cn 




S 




B 


O 


fi 


14H 






0 


a 


o 








CO 


<D 


cx 


CM 


X) 


o 




U 




• 


0) 


cn 


a 


o 


Q) O 


CO 


4J 


C 


B 


o 




rH 




o 


d 


o 






00 


cx 


rH 


o 


CO 


II 




cn 


O 

u 




§ 




O 


0) 



C 

cd 

<D 



a c 

cd 

0 ) Q) 



RMS Vertical Moment = 178152.6 ft-tons 

RMS Peak-to-Peak Vertical Moment = 503,815 ft-tons 

RMS Peak-to-Peak Vertical Bending Stress = 1.99 kpsi 
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The numerical calculation is illustrated below for 
the Roll wind-to-wave relation. 



s 





- ”v)^ ) 



Beaufort No 


h(ft) 


m 

V 


(ton~f t) 


s^^Ocpsi) 


s^Ckpsi) 


1 


7.2 


123079. 


3.465x10^ 


0.573 


.239 


2 


7.6 


158320. 


5.616x10^ 


.878 


.296 


3 


8.0 


186316. 


6.36x10^ 


.0992 


.315 


4 


9.0 


229732. 


5.188x10^ 


.0810 


.285 


5 


10.0 


305014. 


7.572x10^ 


.1184 


.344 


6 


12.4 


443627. 


1.322x10^° 


.2064 


.454 


7 


15.0 


541361. 


2.506x10^° 


. 391.' 


'3 


8 


18.2 


803314. 


3.279x10^® 


.5120 


.716 


9 


22.0 


964491. 


5.530x10^° 


.8637 


.929 


10 


26.0 


1156994. 


7.499x10^° 


1.1712 


1.082 



= S, 






— r 

_ Aii ”l 
~ \z J 



2 

tan e 



Beaufort No 


h(ft) 


s^^Ckpsi) 


H(ft) 


tan 41 ^) 


s^(kpsi) 


(kpsi) 


1 


7.2 


.0573 


.8 


.115692 


2.5 


.374 


2 


7.6 


.0878 


.4 


.115692 


2.5 


.414 


3 


8.0 


.0992 


.6 


.115692 


2.5 


.435 


4 


9.0 


.810 


1.0 


.115692 


2.8 


.430 


5 


10.0 


.1184 


2.5 


.115692 


3.4 


.516 


6 


12.4 


.2064 


2.5 


.115692 


4.3 


.668 


7 


15.0 


.3912 


3.5 


.115692 


5.1 


.758 


8 


18.2 


.5120 


3.8 


.115692 


6.1 


.997 


9 


22.0 


.8637 


4.0 


.115692 


7.0 


1.225 


10 


26.0 


1.1712 


4.0 


.115692 


7.8 


1.403 
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